AD-A094  668  STANFORD  UNIV  CA  DIV  OF  APPLIED  MECHANICS  F/G  11/4 

A  ONE-DIMENSIONAL  EFFECTIVE  DISPERSION  THEORY  FOR  LAYERED  COMPO— ETC(U> 
NOV  80  T  J  DELPH»  G  HERRMANN  N00014-76-C-0054 

UNCLASSIFIED  SUDAM-80-7  NL 


AD  A094868 


»  * 


Contract  N00014-76-C-0054 


A  ONE-DIMENSIONAL  EFFECTIVE  DISPERSION 
THEORY  FOR  LAYERED  COMPOSITES 


by 

T.  J-  Delph 
and 

G.  Herrmann 


SUDAM  Report  No.  80-7 
November  1980 


This  document  has  been  approved  for 
public  release  and  sale;  its  distribution 
if  unlimited. 


o 


81 


■*/<?.!  /  i  .  ?  /  /  / 

UNCLASSIFIED _ _ 


SECURITY  CLASSIFICATION  of  This  page  Daf«  En  farad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

2.  GOVT  ACCESSION  NO. 

0D  -AciiBb? 

3.  RECIPIENT'S  CATALOG  NUMBER 

ii  TITLE  (and  Subtitle)  — *  - 

7^)  b.  9NE-gIMENSI0NAL  EFFECTIVE  DISPERSION  THEORY  fc 
TOR  |AYER£D  COMPOSITES t  '  J 

S.  TYPE  OF  REPORT  ft  PERtO^  COVERED 

’  INTERIM  ,  l  '  j 

s  5:.  r  J  ' 

t.  PERFORMING  ORG.  REPORT  NUMBER 

i  AUTHOR.,) 

^  1  T.  J.IQELPH  aati  G./HERRMANN  \ 

•.  CONTRACT  OR  GRANT  NUMBER^,) 

‘)  n6O914-76-C-(|054  V 

9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

DIVISION  OF  APPLIED  MECHANICS 

STANFORD  UNIVERSITY 

STANFORD,  CALIFORNIA  94305 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  6  WORK  UNIT  NUMBERS 

/  ~  -  / 

NR  064-525  1  — 1 

It.  CONTROLLING  OFFICE  NAME  ANO  AOORESS  '  ,  i 

OFFICE  OF  NAVAL  RESEARCH  J : 

DEPARTMENT  OF  THE  NAVY  V 

ARLINGTON,  VIRGINIA  22217 

-OArTB  / 

pffovairio 0  J 

13.  NUMBER  OF  PAGES 

33 

14.  MONITORING  VGENCY  NAME  S  AOORESSfl/  dlltarent  /root  Cantrolling  Oltice ) 

IS.  SECURITY  CLASS,  (of  thia  report) 

UNCLASSIFIED 

15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thia  Report) 

APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED. 

17.  DISTRIBUTION  STATEMENT  (ot  thm  abstract  entarad  In  Bloc*  20.  II  dlllarent  trom  Raport) 

10.  supplementary  notes 

19.  KEY  WOROS  (Continue  on  raearaa  aide  if  naeaaaary  end  identify  by  block  nvanber) 

COMPOSITES 

LAMINATES 

WAVE  PROPAGATION 

DISPERSION 

2O.^y^0STR ACT  fConf/nti®  on  ravarsa  aide  It  neceaamry  and  /denr//y  by  block  nurr.ber) 

*  A  new  approximate  theory  for  wave  propagation  in  a  periodically-layered 
elastic  body  is  proposed.  The  theory  is  based  on  microstructural  theories  of 
elasticity,  and  makes  use  of  a  collocation  technique  in  conjunction  with  the 
exact  solution  for  wave  propagation  in  an  unbounded  elastic  body.  The  result:; 
are  in  excellent  agreement  with  the  exact  solution.  The  new  theory  is  con¬ 
sidered  to  hold  considerable  promise  for  formulating  approximate  solutions  to 
boundary-initial  value  problems  for  composite  bodies. 

- - - - 1 

1473  coition  of  t  nov  ss  is  obsolete  UNCLASSIFIED' 

/-/  /  '  ■  ■ 


on  Fomu 

WW  |  JAN  73 


S6CURITV  CLASSIFICATION  OF  This  PAGE  r'Ah*n  D*r«  Entartd) 


\ 


ABSTRACT 

A  new  approximate  theory  for  wave  propagation  in  a  periodically-layered 
elastic  body  is  proposed.  The  theory  is  based  on  microstructural  theories  of 
elasticity,  and  makes  use  of  a  collocation  technique  in  conjunction  with  the 
exact  solution  for  wave  propagation  in  an  unbounded  elastic  body.  The  results 
are  in  excellent  agreement  with  the  exact  solution. 

The  new  theory  is  considered  to  hold  considerable  promise  for  formulat¬ 
ing  approximate  solutions  to  boundary-initial  value  problems  for  composite 
bodies. 


1. 


INTRODUCTION 
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The  increasing  use  of  layered  composites  in  light-weight,  high-strength 
structures  has  led  to  a  corresponding  increase  in  efforts  to  analyze  the 
dynamic  behavior  of  these  materials.  In  most  of  these  analyses,  the  com¬ 
posite  is  modeled  as  a  periodically  layered  body,  with  the  constituent 
layers  composed  of  homogeneous,  linearly  elastic  material.  The  layers  may 
be  assumed  to  be  isotropic  or  to  have  some  specified  anisotropy.  Initial/ 
boundary  value  problems  may  then  be  formulated  for  such  a  body  within  the 
framework  of  the  classical  theory  of  elasticity,  but  the  heterogeneous  na¬ 
ture  of  the  body  makes  exact  solutions  very  difficult  to  obtain.  For  this 
reason,  an  approximate  theory  is  usually  require...  for  practical  applications. 

If  the  wavelengths  of  interest  are  considerably  greater  than  the  char¬ 
acteristic  layer  thicknesses,  then  the  effective  modulus  theory  may  be  used. 
Here  the  body  is  modeled  as  a  homogeneous,  but  anisotropic,  solid,  with  the 
material  constants  derived  in  some  fashion  from  the  properties  of  the  con¬ 
stituent  layers.  The  effective  modulus  theory  becomes  invalid,  however, 
when  the  wavelengths  are  on  the  order  of  the  layer  thicknesses.  At  these 
shorter  wavelengths,  the  successive  reflection  and  refraction  of  waves  from 
layer  interfaces  tends  to  break  down  the  structure  of  progressive  waves. 

This  leads  to  a  marked  dispersive  effect  not  predicted  by  the  effective 
modulus  theory. 

To  account  for  this  dispersive  effect,  an  approximate  theory  which  in 
some  fashion  takes  into  account  the  heterogeneous  nature  of  the  body  is  re¬ 
quired.  A  number  of  such  theories,  which  could  be  collectively  called  micro - 
structural  theories,  have  appeared  in  the  literature  in  recent  years.  A 
representative  bibliography  is  given  in  [1] . 

In  this  paper,  we  set  forth  a  new  microstructural  theory  for  wave  prop¬ 


agation  normal  to  the  layering  in  a  periodically  layered  solid.  The  new 


theory,  to  be  known  as  the  effective  dispersion  theory ,  is  conceptually 
simple,  yet  seems  to  be  significantly  more  accurate  than  other  existing 
theories. 

The  effective  dispersion  theory  is  based  in  part  upon  concepts  of 
linear  elasticity  with  microstructure  set  forth  in  the  work  of  Mindlin  [2] 
and  Kunin  [3].  In  the  general  theories  of  these  authors,  no  definite  pro¬ 
cedure  is  given  for  determining  the  coefficients  in  the  differential  equa¬ 
tions  of  motion.  In  the  effective  dispersion  theory,  these  coefficients  are 
determined  by  matching  as  closely  as  possible  the  resulting  approximate 
dispersion  spectrum  to  the  exact  dispersion  spectrum  for  wave  propagation  nor1 
mal  to  the  layering  in  an  unbounded  body.  It  is  from  this  matching  procedure 
that  the  effective  dispersion  theory  takes  its  name.  In  addition,  it  is  pos¬ 
sible  to  match  the  approximate  and  exact  mode  shapes  at  the  end  of  the  first 
Brillouin  zone. 

The  additional  requirements  that  the  approximate  spectrum  have  real 
frequencies  for  all  values  of  wave  numbers ,  and  that  the  coefficients  of  the 
equations  of  motion  be  real,  serve  to  completely  fix  the  coefficients  in  the 
approximate  equations  of  motion.  These  additional  requirements,  however,  im¬ 
pose  certain  restrictions  upon  the  classes  of  materials  for  which  the  effec¬ 
tive  dispersion  theory  is  valid. 

The  resulting  theory  represents  a  significant  improvement  over  most 
existing  approximate  theories.  Not  only  is  it  quantitatively  valid  for 
much  shorter  wavelengths  than  most  of  the  earlier  theories,  but,  perhaps 
more  importantly,  it  embodies  the  first  stopping  band  of  the  spectrum,  i.e. 
the  complex  branch  between  the  first  and  second  Brillouin  zones.  Thus  the 
filtering  property  of  a  periodic  structure  is  contained  in  the  theory,  even 
though  periodicity  as  such  is  no  longer  present  in  the  approximate  spectrum. 
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2.  EQUATIONS  OF  MOTION 

The  layered  solid  is  assumed  to  be  composed  of  periodically  alternat¬ 
ing  layers  of  homogeneous,  isotropic,  linearly  elastic  material  (Fig.  1). 

The  layers  have  thicknesses  2h  ;  2h'  ,  densities  p  ;  p'  ,  and  elastic  con¬ 
stants  (1  ,  p)  ;  (X'  ,  p')  .  Any  two  adjacent  layers  define  a  unit  cell, 
which  is  repeated  periodically  along  the  positive  and  negative  y-axes.  For 
our  purposes,  the  body  will  be  considered  to  be  of  unbounded  extent  along 
the  x-  and  z-axes,  and  bounded  in  the  y-direction  by  two  planes  located  at 
y  ■  yQ  and  y  =  y  having  their  normals  parallel  to  the  y-axis. 

For  wave  propagation  normal  to  the  layering,  the  wave  motion  uncouples 
into  longitudinal  and  shear  wave  motions,  [4] ,  which  may  be  analyzed  sepa¬ 
rately.  In  the  interest  of  definiteness,  we  will  consider  only  shear  wave 
motions,  as  the  case  of  longitudinal  waves  proceeds  in  exactly  analogous 
fashion. 

For  shear  wave  motions,  the  u-displacenent  parallel  to  the  x-axis 
is  assumed  to  be  the  only  non-zero  component  of  displacement.  Then  the 
displacements  in  the  N-th  unit  cell  are  assumed  to  be  given  by 


u(y,t)  =  uQ(y,t) 


y  =0  +  V(y,t> 

N 


yN=° 


u'  (y ,t) 


uQ(y,t) 


(l) 


Here  u  and  u'  are,  respectively,  the  displacement  components  parallel 

to  the  layering  in  the  softer  (unprimed)  layer  and  the  stiffer  (primed) 

layer.  By  the  stiffer  and  softer  layers,  we  mean,  respectively,  the  layers 

with  the  greater  and  lesser  characteristic  shear  wave  speeds.  I  Here  v 

*  1  *N  is 

a  local  coordinate  with  origin  at  the  midplane  and  |y  =0  indicates  the 

N 

global  y-coordinate  is  to  be  evaluated  at  the  layer  midplane. 

Equation  (1)  follows,  to  some  extent,  corresponding  assumptions  made 

in  Mindlin's  microstructural  theory  [2]  and  the  effective  stiffness  theory 

for  composites  [5].  Here,  however,  we  assume  that  the  displacement  field 
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within  the  stiff er  layers  is  constant,  implying  that  all  the  deformation  is 
confined  to  the  softer  layers.  For  the  typical  composite  used  in  structural 
applications,  the  moduli  of  the  stiff er  layers  are  on  the  order  of  50  to 
100  times  greater  than  those  in  the  softer  layers .  Thus  for  such  materials 
this  assumption  can  be  expected  to  agree  well  with  reality  at  the  lower  fre¬ 
quencies  and  wave  numbers  for  which  the  effective  dispersion  theory  is  valid. 
Following  the  notation  of  [5] ,  we  now  define 


el2  =  U0,y/2 


Y21  =  u0,y  -  * 


(2) 


9221  "  3y  (y21}  UQ,yy  '  *,y 


*221  "  ♦.y 


where  e^  is  known  as  the  macro-strain,  y  the  relative  deformation, 

0221  the  gradient  of  relative  deformation,  and  <221  the  local  Reformation 
gradient.  Here  (  )  indicates  partial  differentiation  with  respect  to  y  . 
Again  following  [5] ,  we  postulate  the  following  strain  energy  function  for  the 
layered  solid 

W  =  I  C1S12  +  C2e!2Y21  +  2  C3Y21  +  2  V221  +  Cs‘2210221  +  I  C69221 


2(Cl/4  +  C2  +  c3}  UQ,y  *  (c2/2  +  C3)  °U0,y  +  2  c3*2 


+  4  (c  -  2c  +  c  )  ip 2  +  (c  -  c  )  u  ^  c_u2 

24  5  6  ,y  56  0,yy  ,y  2  6  0,yy 


(3) 
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Equation  (3)  is  similar  in  a  number  of  respects  to  the  strain  energy  func¬ 
tion  postulated  by  Mindlin  [2J ;  Eq.  (3),  however,  includes  some  higher-order 
terms  not  considered  by  Mindlin  and  also  omits  several  terms  included  in 
Mindlin' s  strain  energy  function. 

We  now  set 


Hamilton's  principle  may  be  written  in  the  form 


1  1 
5  J  /  (T  -  W)  dVdt  +  5  /  W  dt  =  0 
t_  v  t  e 


Here  V  is  the  volume  of  the  solid  and  W  represents  the  work  done  by 

e 

external  forces.  Body  forces  are  assumed  to  be  absent.  Following  [2] ,  we 

take  <$W  in  the  form 
e 
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5W 


/  I  (p.  5u  +  Q,  ~  5u^  +  M,  5$) 

^  L  1  0  *1  3y  0  1  'y=y1 

-  (p~5u  +  Q  ~  6u  +  M  i'j)  1  dA 

oo  o  3y  o  o  y=yri  J 


(8) 


where 


parts 


A  is  a  reference  area  measured  in  the  planes  y  **  >  y^  • 

Substitution  of  (5) ,  (6)  and  (8)  into  (7) ,  followed  by  integration  by 
,  yields  the  equations  of  motion  in  the  form 


b,  u„  -  a,  u„  -  a„u„  +  a„t  +  a  ti  =0 

10  1  0,yy  2  0,yyyy  3  ,y  4  ,yyy 


b_ui  -  a„u^  -  a.u  +  a  +  a  '<i>  -0 

2  3  0,y  4  0,yyy  5  6  ,yy 


(9) 


and  the  accompanying  natural  boundary  conditions  at  y  =  yQ  as 


(-a.u„  +  a*  -  a„u  +  a„i^  +  P  )  6u  =  0  , 

1  0,y  3V  2  0,yyy  4  ,yy  0J  0  ' 


(a_u,_  -  a  .If/  +  C>  )  —  cu„  =  0 

2  0,yy  4  ,y  ~Cr  3y  0 


(-a  u  +  a .\l>  +  M.)  t.  =  0  . 

4  0  ,yy  6  ,y  O' 


(10) 


A  similar  set  of  boundary  conditions  applies  at  y  =  y^  . 

Equations  (9)  are  identical  in  form  to  equations  of  motion  for  an  elas¬ 
tic  body  with  microstructure  derived  by  Kunin  [3]  on  the  basis  of  an  analysis 
of  a  generalized  continuum  model.  These  equations  are  considered  by  Kunin  to 
represent  the  most  general  second-order  approximation  to  the  equations  of 
motion  of  an  elastic  body  with  microstructure. 

We  now  wish  to  consider  the  case  of  wave  propagation  in  an  unbounded 
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body.  We  thus  assume  harmonic  solutions  of  the  form 


Ae 


i (ky-wt) 


[>i  =  Be 


i (ky-wt) 


Substitution  into  (9)  leads  to 

(k2a,  -  k4a„  -  ai2b,)  A  +  ik(a,  -  a.k2)  B  =  0 
'1  2  1'  '•34'' 

-  ik(a3  -  a4k2)  A  +  (a^  -  agk2  -  u2b2)  B  =  0  . 


(11) 


(12) 


The  requirement  that  the  determinant  of  the  system  of  Eqs.  (12)  must  vanish 
then  yields  the  dispersion  equation  in  the  form 


aw4  +  (6k4  +  yk2  +6)  w2  +  £k2  +  0k4  +  k6  =  0  , 


(13) 


where 


a  - 


b  b 
12 


S  = 


a_b, 
2  2 


y  ~ 


a6bl  "  alb2 


a,ac  -  a. 
2  6  4 


a2a6  -  a4 


2  6  4 


5  = 


-a5bl 


a  a 
15 


9  = 


2a, a.  -  a,a, 
3  4  16 


a2a5 


(14) 


a,a_  “  a„ 
2  6  4 


a2a6 


a  a  -  a 
2  6  4 


In  both  the  microstructural  theories  of  Minalin  and  Kunin,  no  prescrip¬ 
tion  is  given  for  determining  the  coefficients  a  ,...a_,b  ,b  in  the 

i  biz 

equations  of  motion.  It  is  to  this  question  that  we  address  ourselves  in 
the  succeeding  sections. 


A 
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3.  DISPERSION  SPECTRUM  MATCHING 

In  the  absence  of  analytical  error  estimates,  the  validity  of  an  approx¬ 
imate  solution  in  the  theory  of  elasticity  is  often  judged  by  applying  the 
approximate  theory  to  a  problem  for  which  an  exact  solution  is  known.  The 
degree  of  agreement  between  exact  and  approximate  solutions  is  then  taken  to 
give  a  reasonable  estimate  of  the  error  inherent  in  the  approximate  theory. 


This  idea  of  error  estimation  by  comparison  between  exact  and  approx¬ 
imate  solutions  supplies  the  rationale  for  the  most  important  step  in  the 
construction  of  the  effective  dispersion  theory.  We  will  seek  to  minimize 
the  error  in  the  approximate  solution  by  requiring  the  approximate  disper¬ 
sion  spectrum  to  conform  as  closely  as  possible  to  the  exact  spectrum  for 
wave  propagation  normal  to  the  layering  in  an  unbounded  body.  This  require¬ 
ment  may  be  imposed  in  several  different  ways.  One  way,  for  example,  would 
be  to  minimize  the  mean  square  difference  between  the  two  solutions.  Here, 
however,  we  will  make  use  of  the  simpler  collocation  method,  where  the 
approximate  spectrum  will  be  matched  to  the  exact  spectrum  in  ordinate  and/or 
slope  at  a  number  of  appropriate  points.  Specifically,  we  will  require  that 

k=kB=uB  »  "jk=kB  =  “D 


do 

dk 


u)=0 

k=0 


do 

dk 


ik=kB 


0 


dk 

lk=k3 


0  . 


(15) 


Here  Cg  '  ^3  '  U3  anG  are  in  Fic.  2  and  are  parameters  of  the 

exact  dispersion  spectrum  for  propagation  normal  to  the  layering  in  an  unbounded 
body.  The  equation  defining  the  exact  spectra-,  may  be  found  in  [4] . 
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Upon  introducing  (13)  into  (15),  we  obtain  five  linear  equations  for 
the  six  parameters  ct  ,  6  ,  Y  ,  6  ,  £  ,  and  6  ,  The  five  equations  serve 
to  determine  five  of  the  parameters  in  terms  of  c^  ,  wB  >  u>D  ,  and  k0  , 
and  the  remaining  parameter,  which  we  will  take  to  be  0  •  We  obtain 


a 


2  2 


CJ  M 

B  D 


6  =  (28  +  3kg) 

c 

g 


(16) 


5  -  -  k^  (29  +  3Kg) 


Y  =  -  2 


3k' 


T  +  2Rkn 
2  B 


where 


R  =  - 


(  2  2y 

2  2 
B  D 
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4.  DETERMINATION  OF  9 

In  order  to  completely  determine  the  approximate  dispersion  spectrum 
(13) ,  we  must  specify  a  value  of  9  .  To  this  end,  we  will  first  require 
that  the  approximate  spectrum  have  real-valued  frequencies  for  all  values  of 
wave  number.  This  requirement  insures  the  stability  of  the  displacements 
(1) ,  in  that  exponential  growth  of  the  displacements  with  time  is  not 
allowed.  From  (13) , 

w2  =  _  (8k4  +  yk2  +  6)  ±  \(gk4  +  yk2  +  6 ) 2  -  4ak2(€  +  9k2  +  k4) 

W  2a 


From  (17) ,  necessary  and  sufficient  conditions  that  u>  be  real-valued  are 


(8k4  +  yk2  +  6)  2  -  4ctk2(£  +  6k2  +  k4)  >  0 


(18) 


and 


ct(£  + 


0k2  +  k4) 


^  0 


(8k4  + 


2 

yk  + 


6) 


2a 


>  0 


(19) 


Examining  the  first  of  (19) ,  substitution  from  (16)  results  in 


(0  +  2kf) [k4  +  0k2  -  k2  (29  +  3k2)]  <  0  . 
B  B  B 


(20) 


Thus  the  two  factors  in  (20)  must  be  of  opposite  sign  or  vanish  for  all 
values  of  k  .  If  we  denote  the  second  factor  by 


g(k) 


4  2  2  2 

=  k  +  rk  -  k^(2r  +  3k„)  , 

B  B 


(21) 


then  the  roots  of  g(k)  are  given  by 
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t2  -  -  ~  ±  W(9  +  2k2)  (8  +  6k2)  . 


(22) 


Now  g(k)  will  maintain  the  same  sign  for  all  values  of  k  ,  or  at  most 
vanish  at  an  isolated  value  of  k  ,  if  it  has  only  complex  roots,  or  at  most 
a  double  real  root.  From  (22) ,  this  will  be  the  case  if 


-6kg  <  9  <  - 


2k_ 


(23) 


2 

However,  from  (16) ,  9  =  -2k  implies  a  =  0  ,  which  results  in  a  degener- 

B 

ate  dispersion  equation.  Thus  we  reject  this  value  and  rewrite  (23)  as 


2  2 

-6k  i  9  <  -2k^ 
D  B 


(24) 


Then  for  values  of  9  satisfying  (24),  g(k)  can  be  seen  to  be  positive,  or 
zero  for  all  real  values  of  k  ,  while  the  first  factor  in  (20)  is  negative, 
thus  satisfying  the  inequality. 

We  now  consider  the  second  of  (19) .  Since  a  >  0  for  9  in  the  range 
(24) ,  the  second  of  (19)  implies  that 


8k4  +  yk2  +  5  <  0  . 


(25) 


Now  o  <  0  if  (24)  is  satisfied,  so  the  inequality  (25)  will  be  satisfied 

for  all  values  of  k  if  and  only  if  the  biquadratic  expression 
4  2 

8k  +  yk  +5=0  has  no  real  roots,  or  at  most  a  double  real  root.  This 
requirement  implies  that  the  discriminant  must  satisfy  the  inequality 


Y2  -  485  s  0  . 


(26) 


If  we  set 


,  where 


2  <  c  S  6  (27) 

from  (24)  and  substitute  from  (16) ,  (26)  becomes 

f(c,Rk2c2)  i  (2  -  c)  Rk2c2  +  (3  -  2c)  <  0  (28) 

2  2 

where  we  note  that  Rk  c  <  0  .  The  region  of  common  validity  of  the  in- 

b  g 

equalities  (28)  and  (27)  is  shown  as  the  hatched  area  in  Fig.  3.  It  is 

apparent  from  the  figure  that  we  can  always  pick  a  value  of  c  so  as  to 

2  2  9 

satisfy  both  (27)  and  (28) .  We  note  in  addition  that  for  Rc  k  <  -  —  ,  (28) 

g  B  4 

represents  a  more  severe  restriction  on  the  upper  limit  of  c  than  does  (27) 
We  now  examine  Eq.  (18),  which  may  be  written  as 

h (k)  =  S2k8  +  (20y  -  4a)  k6  +  (y2  +  266  -  4a0)  k4 

+  (2y6  -  4a£)  k2  +  S2  k  0  (29) 

2 

h(k)  is  a  fourth-order  polynomial  in  k  ,  whose  coefficients  are  functions 
of  the  parameter  9  through  (16).  Since  h(0)  =6  >  0  ,  it  is  apparent 
that  a  necessary  and  sufficient  condition  that  (29)  hold  for  all  values  of 
k  is  that  h(k)  have  no  real  roots,  or  at  most  double  real  roots. 

Due  to  the  complexity  of  h(k)  ,  numerical  ar.alysis  is  required  to 
determine  whether  or  not  there  exist  values  of  9  such  that  h(k)  has  no 
real  roots,  or  double  real  roots,  and  such  that  (24)  and  (28)  are  satisfied. 
If  this  is  the  case,  then  for  these  values  of  9  ,  the  dispersion  equation 


(13)  will  yield  only  real  frequencies  and  thus  guarantee  stable  displacements 
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To  complete  the  determination  of  the  parameter  0  ,  we  now  consider 
the  question  of  inverting  equations  (14)  to  determine  the  coefficients  of 
the  equations  of  motion  in  terms  of  the  six  parameters  a  ,  8  ,  y  *  $  *  £  # 
and  9  .  Since  the  equations  of  motion  are  homogeneous,  one  of  the  coef¬ 

ficients  may  be  fixed  arbitrarily.  We  will  take  this  coefficient  to  be  a^ 
and  normalize  the  remaining  seven  coefficients  with  respect  to  a^  .  The 
inversion  procedure  then  yields 


where  the  ratio  a  /a, 

J  o 


is  a  root  of 
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The  positive  sign  in  the  last  term  in  (31)  is  to  be  taken  when  and  a 

are  of  the  same  sign,  and  the  negative  sign  when  they  are  of  the  opposite 

sign.  The  ratio  a./a,.  is  left  as  yet  undetermined. 

2  o 

Since  9  is  not  completely  determined  by  the  requirement  that  the  fre¬ 
quencies  in  the  approximate  spectrum  be  real,  we  now  seek  to  obtain  a  value 
of  9  which,  besides  yielding  real  frequencies,  also  gives  real  roots  for 
a^/a^  in  (31) .  Although  there  seems  to  be  no  a  priori  requirement  that 
a^/a^  ,  or  any  of  the  other  coefficients  in  the  equations  of  motion,  be  real, 
the  imposition  of  this  requirement  will  simplify  the  resulting  theory.  Also, 
it  will  permit  the  determination  of  a  unique  value  of  9  ,  as  will  shortly 
be  demonstrated. 

Squaring  (31) ,  we  obtain  the  following  fourth-order  polynomial  in  a_/a, 

5  6 


whore  me  coefficients  are  functions  of  the  parameter  f  .  Due  to  its  com- 
,  tin  root:,  of  id  r.usr.  hi-  analyze!  numerical  ip. 


When  this 


done , 
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a  rather  surprising  fact  emerges,  namely,  that  (32)  and  h(k)  as  defined 
by  (29)  both  exhibit  a  double  real  root  for  the  same  value  of  0  .  Further¬ 
more,  for  values  of  0  greater  than  this  value,  (32)  and  h(k)  both  have 
distinct  real  roots,  while  for  values  less  than  this  value,  they  have  no 
real  roots. 

The  condition  on  the  coefficients  of  a  polynomial  such  that  the  poly¬ 
nomial  has  double  real  roots  may  be  determined  from  the  theory  of  eliminants 

4  3  2 

[6] .  For  a  fourth-order  polynomial  of  the  form  aQx  +  a^x  +  a2x  +  a3x 
+  a4  =  0  ,  the  criterion  for  double  roots  may  be  expressed  as 

(3ai  '  8a0a2)  (2ai32  ‘  12a0a3)  (aia3  '  l6a0a4) 

(4a2  -  8a1a;J  -  16aQa4)  <2a2a3  “  12a1a4)  =  0 

Ula3  “  l6a0a4}  {2a2a3  ‘  l2aia45  (3a3  '  ®a2a4) 

Equation  (33)  ,  when  applied  to  (32)  and  h(k)  ,  becomes  too  complex  to  be 

treated  analytically,  and  must  be  evaluated  numerically.  This  was  done  for 

several  different  sets  of  values  of  u>„  ,  ,  c  ,  and  k  .  In  every  case 

B  D  g  B 

the  same  value  of  9  produced  a  double  real  root  for  both  h(k)  and  (32) . 

It  is  thus  concluded  that  this  behavior  is  not  coincidental. 

We  now  see  that  0  may  be  determined  uniquely  as  that  value  which 
yields  double  real  roots  for  h(k)  and  (32)  ,  thus  satisfying  simultaneously 
the  inequality  (29)  and  the  condition  that  a  /a.  be  real.  Furthermore, 

D  O 

the  quantity  a./a^  will  be  determined  uniquely  as  well.  The  ambiguity 
in  the  signs  of  a^  and  a4  may  be  resolved  by  numerically  determining 
whether  the  positive  or  negative  sign  must  be  taken  in  the  last  term  of  (31) 
in  order  that  the  double  real  root  a  /a  satisfy  (31).  Of  course,  the  remaining 

J  O 

two  roots  of  h(k)  must  be  complex  to  insure  real  frequencies.  Numerical 
analysis  is  required  to  test  the  satisfaction  of  this  requirement. 


Once  0  is  determined,  the  approximate  dispersion  spectrum  is  completely 
specified.  As  a  numerical  example,  we  take  u/u '  -  0.02  ,  up'/u'p  =  0.06 
and  h'/h  *  4  ,  and  introduce  a  nondimensional  frequency  and  wave  number  de¬ 
fined  by  Q  *  2hu)/(ir*/u/p)  and  n  =  (2h/n)k  .  For  these  values,  the  exact 

solution  [4]  yields  Q_  =  2hu)D/(ir/U7p)  =  0.179005  ,  Q  =  0.831029  , 

o  B  D 

c  H  (dSl/dn)  I  n  =  1.334401  ,  and  =  (2h/Tr)k_  =  0.20  .  We  now  apply 
y  j  n=U  B  B 

the  matching  procedure  with  these  nondimensional  parameters.  The  value  of 

9  which  yields  double  real  roots  for  h(k)  and  (32)  is  found  to  be 

2 

9  -  -0.156317  =  -3.90793n„  ,  which  satisfies  (24).  The  remaining  pair  of 

B 

2  2 

roots  of  h(k)  is  found  to  be  a  complex  conjugate  pair.  Also,  Rc^nB  =  -2.32594, 
satisfying  (28) . 

The  resulting  spectrum  is  compared  to  the  exact  spectrum  for  shear 
waves  propagating  normal  to  the  layering  in  Figs.  4  and  5.  Figure  4  shows 
the  exact  and  approximate  spectra  over  the  first  two  Brillouin  zones  for 
real  values  of  the  wave  number  n  ,  plotted  on  the  periodic  zone  scheme. 

Since  only  waves  with  non-negative  group  velocity  are  to  be  considered,  only 
the  acoustic  branch  in  the  first  Brillouin  zone  and  the  optical  branch  in 
the  second  Brillouin  zone  are  of  interest.  It  can  be  seen  that  the  approxi¬ 
mate  spectrum  agrees  extremely  well  with  the  exact  solution.  The  spectrum 
of  a  typical  existing  approximate  theory,  the  effective  stiffness  theory  [5] , 
is  shown  for  comparison.  It  is  seen  that  the  effective  dispersion  theory 
yields  a  substantial  improvement  in  accuracy.  Figure  5  shows  the  exact  and 
approximate  complex  branches  connecting  the  acoustic  and  optical  branches. 

Here  again  the  agreement  between  the  two  solutions  is  quite  good. 

In  passing,  it  should  be  noted  that  the  requirement  that  h(k)  have  a 
double  real  root  implies  that  the  optical  and  acoustic  branches  of  the 
approximate  spectrum  will  touch  at  a  point.  For  the  parameters  used  in  the 
numerical  example,  this  point  lies  in  the  third  Brillouin  zone. 
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5.  MATCHING  OF  MODE  SHAPES 

We  now  seek  to  determine  the  quantity  a  /a  by  imposing  the  require- 

4  O 

ment  that  the  approximate  mode  shape  match  the  exact  mode  shape  as  closely 
as  possible  at  the  right-hand  end  of  the  first  Brillouin  zone. 

In  order  to  compute  the  approximate  mode  shapes,  we  arbitrarily  assume 
the  global  coordinate  system  to  have  its  origin  at  the  midplane  of  the  layer 
with  primed  constants  in  the  zeroth  unit  cell.  Using  this  scheme,  the  mid¬ 
plane  of  the  layer  with  primed  constants  within  the  N-th  unit  cell  is  located 
at 


y  =  2N(h  +  h') 


(34a) 


while  the  midplane  of  the  layer  with  unprimed  constants  is  located  at 


y  =  (1  +  2N) (h  +  h')  .  (34b) 


Thus,  introducing  (11)  into  (1)  and  making  use  of  the  nondimens ional  vari¬ 
ables  defined  in  the  previous  section,  the  approximate  displacement  mode 
shapes  (1)  take  the  form 


u 

A 


1 


i  (2N+1)  (1+e) "h/2 
e 


u' 

A 


iN (l+:)~r 
e 


(35) 


£_noe  we  are  usino  complex  :iotac_nr. .  it  is  implicitly  understood  that  only 
tne  real  oart  if  ;.r  wt—i  oe  ::f  s_  or.  if  trance .  Here  the  time  variation  has 


- —4 

4 
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been  omitted,  and  N  *  0  ,  ±1  ,  ±2  ,  . . .  .  It  may  be  seen  immediately  that 

the  displacements  defined  by  (35)  have  the  same  periodicity  as  the  exact 

2hkB  1 

solution  at  the  end  of  the  first  Brillouin  zone  where  ng  ■  — ~ —  =  > 

i.e. , periodic  with  a  period  of  two  unit  cells. 

The  ratio  (Bh/A)  may  be  determined  from  the  first  of  (12)  ,  which  gives 


Bh 

A 


2  4 

n  al  “  n  a2  “ 
3 

na,  -  n 


(36) 


where  the  constants  involved  in  the  nondimensionalization  procedure  are 

assumed  to  have  been  absorbed  into  the  parameters  a^  ,  . . .  ,  b^  .  We  now 

write  (Bh/A)  *  ic  ,  and  evaluate  the  displacements  (35)  at  r)  =  n  • 

B 

We  obtain 


(37) 


For  the  parameters  used  in  the  numerical  example  in  the  previous  section, 
the  displacement  corresponding  to  the  exact  solution  is  symmetric  about  the 
layer  midplanes  in  the  layers  with  primed  constants  and  antisymmetric  in 
the  layers  with  unprimed  constants.  This  pattern  of  symmetry  can  be  seen 
to  be  realized  in  the  approximate  displacements  (37) .  If  we  now  require 
continuity  of  displacement  at  the  interface  of  the  N-th  unit  cell,  where 
y^/h  =  -1  ,  we  obtain 
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C  =  1  . 


(38) 


Using  (36) ,  (38) ,  and  (30) ,  we  may  solve  for  the  required  value  of  a2/afe 
which  is 


(39) 


The  quantity  a_/a.  is  determined  in  the  previous  section.  The  positive 

D  O 

sign  is  taken  in  the  second  term  of  the  numerator  of  (40)  if  a,  and  a. 

3  4 

are  taken  to  be  of  opposite  sign,  while  the  negative  sign  is  taken  if 
and  a4  are  of  the  same  sign. 

The  exact  and  approximate  displacement  mode  shapes  at  the  right-hand 

end  of  the  first  Brillouin  zone  are  sketched  in  Fig.  6  for  a  /a  defined 

Z  o 

by  (39).  The  exact  mode  shapes  are  taken  from  [7] .  Here  l  and  Z ’  are 
local  coordinates  having  their  origins  at  the  midplanes  of  the  layers  with 
unprimed  and  primed  constants,  respectively,  and  being  normalized  by  the 
layer  thickness.  The  agreement  between  exact  and  approximate  mode  shapes 
in  Fig.  6  can  be  seen  to  be  excellent.  Figure  7  shows  the  exact  and  approx¬ 
imate  mode  shapes  at  two  points  in  the  interior  of  the  first  Brillouin  zone. 
Here  the  agreement  is  not  as  good  as  at  the  zone  end.  This  is  not  surprising 
in  view  of  the  fact  that  the  exact  displacements  are  in  general  quasi- 
periodic,  as  discussed  in  [7] ,  while  the  approximate  displacements  given 
by  (1)  and  (11)  are  strictly  periodic.  Also,  continuity  of  displacement  at 
the  layer  interfaces  in  the  effective  dispersion  theory  is  only  approximate. 


as  may  be  noted  from  the  figure. 
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6.  DISCUSSION  AMD  CONCLUSIONS 

As  mentioned  earlier,  there  exist  in  the  literature  a  number  of 
approximate  theories  for  wave  propagation  in  a  layered  composite  of  the 
type  considered  in  this  work.  For  the  case  of  wave  propagation  normal 
to  the  layering,  these  theories  are  generally  limited  in  validity  to 
regions  in  the  interior  of  the  first  Brillouin  zone.  The  ability  to 
model  this  region  of  the  spectrum  accurately  is  certainly  the  most  im¬ 
portant  consideration  for  stress  analysis  ;  however,  it  seems  likely 
that  the  stopping  band  between  the  first  and  second  Brillouin  zones  may 
also  play  a  significant  role  in  some  practical  problems.  The  effective 
dispersion  theory  outlined  in  this  paper  is  the  only  existing  approxi¬ 
mate  theory  which  has  been  shown  to  contain  this  stopping  band.  In  ad¬ 
dition,  the  effective  dispersion  theory  is  more  accurate  than  most  other 
proximate  theories  over  the  first  Brillouin  zone.  It  is  quite  accurate 
as  well  over  the  second  Brillouin  zone,  a  region  which  other  approximate 
theories  have  not  attempted  to  model.  For  these  reasons,  the  effective 
dispersion  theory  is  thought  to  represent  a  significant  advance  in  the 
construction  of  approximate  theories  for  wave  propagation  in  layered 
solids. 

The  effective  dispersion  theory  is  subject  to  one  limitation,  how¬ 
ever,  in  that  there  seems  to  be  no  guarantee  that  a  value  of  0  can 
be  found  so  that  inequalities  (24) ,  (23)  ,  and  (29)  are  simultaneously 
satisfied.  It  appears  that  these  inequalities  must  be  individually 
tested  for  each  set  of  material  and  geometric  parameters  under  consider¬ 
ation,  and  that  no  definite  general  statement  can  be  made  as  to  the 
range  of  applicability  of  the  effective  dispersion  theory. 

On  the  other  hand,  the  numerical  example  presented  in  this  paper 
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was  chosen  to  represent  a  typical  composite  material,  and  here  it  was 
shown  that  it  was  possible  to  satisfy  all  the  relevant  inequalities. 

It  therefore  seems  likely  that  the  effective  dispersion  theory  will  be 
applicable  to  the  analysis  of  a  number  of  classes  of  composite  materials 
of  practical  interest. 

In  closing,  we  will  briefly  discuss  two  areas  where  further  effort 
is  required  to  refine  the  effective  dispersion  theory.  The  first  of 
these  is  the  question  of  uniqueness  of  solution.  It  has  been  shown  in 
[5]  that  uniqueness  of  solution  follows  from  the  positive  definiteness 
of  the  strain  and  kinetic  energy  densities.  More  recently,  Edelstein 
[8]  has  proved  uniqueness  of  solution  for  microstructural  theories  of 
the  Mindlin  type  if  the  strain  energy  density  is  homogeneous  of  order 
two  and  if  the  coefficients  in  the  kinetic  energy  density  expression, 
and  b2  ,  sire  of  the  same  sign.  Unfortunately,  because  the  effec¬ 
tive  dispersion  theory  places  primary  emphasis  upon  matching  procedures, 
it  is  not  possible  to  satisfy  these  conditions.  Thus,  the  question  of 
uniqueness  of  solution  remains  an  open  one. 

One  rather  general  criterion  for  uniqueness  may  be  stated  immedi¬ 
ately  however.  It  has  been  shown  [9]  that  stability  in  the  sense  of 
Liapunov  represents  a  necessary  and  sufficient  condition  for  unique¬ 
ness  of  solution  in  dynamic  systems.  Thus,  if  a  Liapunov-stable  solu¬ 
tion  can  be  found  to  equations  (9) ,  subject  to  the  boundary  conditions 
(10)  and  an  appropriate  set  of  initial  conditions,  the  uniqueness  of 
the  solution  is  insured. 

The  second  area  requiring  further  effort  concerns  the  boundary 
conditions  (10) .  At  the  present  time  we  are  unable  to  give  a  prescrip¬ 
tion  for  calculating  the  boundary  forces  P  ,  Q  ,  and  M  for  a  pure 
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traction  boundary  value  problem.  This  difficulty  seems  to  be  common  to 
theories  of  the  microstructural  type  and  represents  an  impediment  to  the 
solution  of  realistic  boundary-initial  value  problems .  The  most  obvious 
method  for  treating  a  traction  boundary  value  problem,  setting  P  equal  to 
the  applied  traction  and  Q  *  M  =  0  ,  appears  to  be  somewhat  questionable 
here,  since  setting  Q  ■  M  =  0  in  equations  (10)  leads  to  uo,yy  =  ^,y  =  0 
at  the  bounding  surface.  These  conditions  seem  to  be  unnecessarily  restric¬ 
tive.  It  is,  of  course,  possible  to  formulate  and  solve  boundary  value 
problems  of  the  mixed  variety,  but  these  are  seldom  of  practical  interest. 

The  resolution  of  this  question  must  at  the  present  time  await  further  study. 
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FIGURE  CAPTIONS 

Fig.  1  Geometry  of  layered  elastic  solid 

Fig.  2  Parameters  of  exact  dispersion  spectrum  used  in  effective  dis¬ 

persion  theory. 

Fig.  3  Region  of  common  validity  of  inequalities  (27)  and  (28) . 

Fig.  4  Spectrum  of  effective  dispersion  theory  over  the  first  two  Bril- 
louin  zones  (real  branches) . 

Fig.  5  Complex  branch  between  the  first  two  Brillouin  zones. 

Fig.  6  Effective  dispersion  mode  shapes  and  exact  mode  shapes  at  the 

right-hand  end  of  the  first  Brillouin  zone. 

Fig.  7  Effective  dispersion  mode  shapes  and  exact  mode  shapes  at  two 
points  in  the  interior  of  the  first  Brillouin  zone. 
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